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At an hadronic scale the effect of CP- violating interactions that typically appear in extensions of 
the Standard Model may be described by an effective Lagrangian, in which the operators are ex- 
pressed in terms of lepton and partonic gluon and quark fields, and ordered by their mass dimension, 
k > 4. Using a global-symmetry-preserving truncation of QCD's Dyson-Schwinger equations, we 
compute the p-meson's electric dipole moment (EDM), d p , as generated by the leading dimension- 
four and -five CP-violating operators and an example of a dimension-six four-quark operator. The 
two dimension-five operators; viz., quark-EDM and -chromo-EDM, produce contributions to d p 
whose coefficients are of the same sign and within a factor of two in magnitude. Moreover, should 
a suppression mechanism be verified for the S-term in any beyond-Standard-Model theory, the con- 
tribution from a four-quark operator can match the quark-EDM and -chromo-EDM in importance. 
This study serves as a prototype for the more challenging task of computing the neutron's EDM. 

PACS numbers: 11.30.Er, 14.40.Be, 11.15.Tk, 12.60.Jv 



I. INTRODUCTION 

The action for any local quantum field theory is invari- 
ant under the transformation generated by the antiuni- 
tary operator CPT ', which is the product of the inver- 
sions: C, charge conjugation; P, parity transformation; 
and T, time reversal. The combined CPT transforma- 
tion provides a rigorous correspondence between parti- 
cles and antiparticles, and it relates the S matrix for any 
given process to its inverse, where all spins are flipped 
and the particles replaced by their antiparticles. Lorentz 
and CPT symmetry together have many consequences, 
amongst them, that the mass and total width of any par- 
ticle are identical to those of its antiparticle. 

It is within this context that the search for the in- 
trinsic electric dipole moment (EDM) of an elementary 
or composite but fundamental particle has held the fas- 
cination of physicists for over sixty years Its exis- 
tence indicates the simultaneous violation of parity- and 
time-reversal-invariance in the theory that describes the 
particle's structure and interactions; and the violation of 
P- and T-invariance entails that CP symmetry is also 
broken. This last is critical for our existence because 
we represent a macroscopic excess of matter over anti- 
matter. As first observed by Sakharov Q, in order for 
a theory to explain an excess of baryon matter, it must 
include processes that change baryon number, and break 
C- and CP-symmetries; and the relevant processes must 
have taken place out of equilibrium, otherwise they would 
merely have balanced matter and antimatter. (Alter- 
nately, the presence of CPT violation can circumvent the 
out-of-equilibrium environment.) 



The electroweak component of the Standard Model 
(SM) is capable of satisfying Sakharov's conditions, ow- 
ing to the existence of a complex phase in the 3 x 3-CKM 
matrix which enables processes that mix all three quark 
generations. However, this high-order process is too 
weak to explain the observed matter-antimatter asym- 
metry 0-[5|. Hence, it is widely expected that any de- 
scription of baryogenesis will require new sources of CP 
violation beyond the SM. This presents little difficulty, 
however, because extensions of the SM typically pos- 
sess CP-violating interactions, whose parameters must, 
in fact, be tuned to small values in order to avoid con- 
flict with known bounds on the size of such EDMs . 
(For recent analyses, see, e.g., Refs. (Iol-[T2| and references 
therein.) 

The question here is how such bounds should be im- 
posed. That is not a problem for elementary particles, 
like the electron. However, it is a challenge when the 
SM extension produces an operator involving current- 
quarks and/or gluons. In that case the CP violation is 
expressed as an hadronic property and one must have at 
hand a nonperturbativc method with which to compute 
the impact of CP- violating features of partonic quarks 
and gluons on the hadronic composite. 

To elucidate, extensions of the SM are typically active 
at some large but unspecified energy-scale, A, and their 
effect at an hadronic scale is expressed in a low-energy 
effective Lagrangian: 

/•,:- » A ; \ (1) 
where are composite CP-odd local operators of di- 
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mcnsion k > 4 and {Kj} are dimensionless strength pa- 
rameters, which monitor the size of the model's CP- 
violating phases and commonly evolve logarithmically 
with the energy scale. The calculation of an hadronic 
EDM therefore proceeds in two steps. The first, eas- 
ier, part requires calculation of the coefficients {Ki\ in 
a given model. This involves the systematic elimination 
of degrees-of-freedom that are irrelevant at energy-scales 
less than A. The second, far more challenging exercise, is 
the nonperturbative problem of translating the currcnt- 
quark-level interaction in Eq. (JTJ) into observable proper- 
ties of hadrons. 

We illustrate the procedure in the case of the p-meson. 
Not that there is any hope of measuring a p-meson EDM 
but because the nonperturbative methods necessary can 
most readily be illustrated in the case of systems defined 
by two valence-quark degrees-of-freedom. In taking this 
path, we follow other authors fl3l - fl5| but will nonethe- 
less expose novel insights, especially because we consider 
more operator structures than have previously been con- 
sidered within a single unifying framework. It is worth 
remarking here that particles with spin also possess a 
magnetic dipolc moment. That moment is aligned with 
the particle's spin because it is the only vector available. 
The same is true of the expectation value of any electric 
dipolc moment. 

Herein we shall estimate the contribution of some di- 
mension four, five and six operators to the EDM of the 
p + -meson; viz., the impact on the p of the local La- 
grangian density 
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where: latin superscripts represent colour; g s is the 
strong coupling constant; F pu and G^ v arc photon and 

gluon field-strength tensors, respectively, and G a ^ v = 
{l/2)e pvXp G a Xp - {Q t \i = 1,2} = {u L ,d L }, with the sub- 
script indicating left-handed; 9 is QCD's effective 6- 
parameter, which combines #qcd and the unknown phase 
of the current-quark-mass matrix; and {d q }, {d q } are 
quark EDMs and chromo-EDMs, respectively. 

We note that Eq. ((2]) is expressed at a renormalisa- 
tion scale ( ~ 2 GeV, which is far below that of elec- 
troweak symmetry breaking but still within the domain 
upon which pcrturbativc QCD is applicable. Moreover, 
we have chosen to include just one dimension-six oper- 
ator in the Lagrangian; i.e., a particular type of four- 
fcrmion interaction. There is a host of dimension-six op- 
erators, Weinberg's CP-odd three-gluon vertex amongst 
them [16[. However, for our illustrative purpose, nothing 
is lost by omitting them because the potency of the one 
operator we do consider can serve as an indication of the 



strength with which each might contribute. 

One merit of our analysis of the contribution from 
Eq. to the EDM of the p + -meson is the connection 
of these EDM responses with values of a vast array of 
hadron observables that are all computed within precisely 
the same framework using exactly the same parameters 
[1714221 ]. We explain this framework in Sec.HH In addition 
to providing the first such comprehensive treatment, our 
study is novel in considering the impact of a dimension- 
six operator on the p + -meson's EDM. 

We introduce the p-meson electromagnetic form factors 
in Sec. IIII1 The effects of Eq. on the p-meson bound- 
state are analysed in Sec. lIVl Each interaction term is 
considered separately, so that we present a raft of alge- 
braic formulae that are readily combined, evaluated and 
interpreted. Numerical results are provided in Sec.lVland 
placed in context with previous studies. Section [VTl is an 
epilogue. 

II. p-MESON AS A BOUND STATE 

A. p-7 Vertex 

The p + -meson is a composite particle and thus its 
EDM appears in the dressed vertex that describes its 
coupling with the photon; viz., 

= <2 T aa , (p) {(p + p')^-S a , p ,'E(q 2 ) + q a , q p , d(q 2 )] 

-(Sfia'Qp 1 - 5^q al )M{q 2 ) 

-ie a ^'„aqa'D(q 2 )}^p(p'), (3) 

where: p a is the momentum of the incoming p-mcson; 
p'p , that of the outgoing p; q p = p' — p^ ; and 

a5,(p) = 4*-^. (4) 

The vertex involves four scalar form factors whose q 2 = 
values are understood as follows: £(0), electric charge, 
which is "1" in this case; fW(0), magnetic moment, p p . in 
units of e/[2m p ], where e is the magnitude of the electron 
charge; Qjfi) = {2 / m 2 )(Q p + — 1), with Q p the meson's 
electric quadrupolc moment; and 2)(0) is the meson's 
electric dipole moment, in units of e/[2m p ]. 

B. Contact Interaction 

Our goal is calculation of the last of these, ©(0), 
and for this we choose to work within the continuum 
framework provided by QCD's Dyson-Schwinger equa- 
tions (DSEs) (2314251 I . To be specific, we perform the 
computation using a global-symmetry-preserving treat- 
ment of a vector x vector contact-interaction because that 
has proven to be a reliable explanatory and predictive 
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tool for hadron properties measured with probe momenta 
less-than the dressed-quark mass, M ~ 0.4 GeV [T3-H3- 

To expand upon the reasons for this choice of interac- 
tion we note that DSE kernels with a closer connection to 
perturbative QCD; namely, which preserve QCD's one- 
loop renormalisation group behaviour, have long been 
employed in studies of the spectrum and interactions 
of mesons [26T - |28j . Such kernels arc developed in the 
rainbow-ladder approximation, which is the leading-order 
in a systematic and global-symmetry-preserving trunca- 
tion scheme [29l . |30| ; and their model input is expressed 
via a statement about the nature of the gap equation's 
kernel at infrared momenta. With a single parameter 
that expresses a confinement length-scale or strength 
[3TI HH , they have successfully described and predicted 
numerous properties of vector [32H36I ] and pseudoscalar 
mesons [32|, [3a - l40l | with masses less than 1 GeV, and 
ground-state baryons j4lTj44l |. Such kernels are also reli- 
able for ground-state heavy- heavy mesons [45j]. Given 
that contact-interaction results for low-energy observ- 
ables are indistinguishable from those produced by the 
most sophisticated interactions, it is sensible to capitalise 
on the simplicity of the contact-interaction herein. 

The starting point for our study is the dressed-quark 
propagator, which is obtained from the gap equation: 
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9 D^(p- q)-J^S(q) —T v [q,p), (5) 



wherein m is the Lagrangian current-quark mass, D^ v is 
the vector-boson propagator and T u is the quark-vector- 
boson vertex. We use 



2 D^(p-q)=5, 



4ira m 



(6) 



where ma = 0.8 GeV is a gluon mass-scale typical of the 
one-loop renormalisation-group-improved interaction in- 
troduced in Ref. [H[, and the fitted parameter am /ir = 
0.93 is commensurate with contemporary estimates of 
the zero-momentum value of a running-coupling in QCD 
[46l [47} . Equation ^ is embedded in a rainbow-ladder 
truncation of the DSEs, which is the leading-order in 
the most widely used, symmetry-preserving truncation 
scheme [301 . This means 



1^0,9) = lv 



(7) 



in Eq. ([5]) and in the subsequent construction of the 
Bcthc-Salpeter kernels. One may view the interaction 
in Eq. ((6]) as being inspired by models of the Nambu- 
Jona-Lasinio (NJL) type (48|. However, in implementing 
the interaction as an element in a rainbow-ladder trun- 
cation of the DSEs, our treatment is atypical; e.g., we 
have a single, unique coupling parameter, whereas com- 
mon applications of the NJL model have different, tun- 
able strength parameters for each collection of operators 
that mix under symmetry transformations. 



Using Eqs. ((6]), (0, the gap equation becomes 

a-lf \ ■ 1 1 167raiR f d 4 q 

S (p) = 27 • p + m +—-^T J 7m S{q) 7 P , (8) 

an equation in which the integral possesses a quadratic 
divergence, even in the chiral limit. When the divergence 
is regularised in a Poincare covariant manner, the solu- 
tion is 

Sip)' 1 =ij-p + M, (9) 
where M is momentum-independent and determined by 
4a m r. 1 



M = m + M 



3nm G j 



ds s ■ 



M 2 



(10) 



Our regularisation procedure follows Ref. [431; i.e., we 
write 



1 



M 2 







dTe -r ( s+M-) 

" dre-< s+M ^ (11 
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M 2 



(12) 



where - 



are, respectively, infrared and ultraviolet reg- 



ulators. It is apparent from Eq. (|12p that Tj r ~: 1/Aj r fi- 
nite implements confinement by ensuring the absence of 
quark production thresholds [23|, • Since Eq. © docs 
not define a rcnormalisablc theory, then A uv := 1/t uv 
cannot be removed but instead plays a dynamical role, 
setting the scale of all dimensioned quantities. 
Using Eq. (TiT]) , the gap equation becomes 



M = m 



M^^- C iu (M 2 ) . 



(13) 



where 



C in {M 2 ) = M 2 C (M 2 ) 



(14) 

= M 2 [r(-l, A/ 2 r u 2 v ) - T(-l, M 2 t 2 )} , (15) 

with r(a,y) the incomplete gamma-function, and, for 
later use, C[ u (z) = -~z(d/dz)C iu (z). 

In rainbow-ladder truncation, with the interaction in 
Eq. (O, the homogeneous Bethe-Salpeter equation for the 
colour-singlet p-meson is 



* 3 m 



167TOTR 



W) 4 



(16) 



where x£(g;P) = S(q + P)TP^q;P)S(q) and T^P) 
is the meson's Bethe-Salpeter amplitude. Since the 
integrand does not depend on the external relative- 
momentum, k, then a global-symmetry-prescrving reg- 
ularisation of Eq. (|16[) yields solutions that are indepen- 
dent of k. With a dependence on the relative momentum 
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TABLE I. Results obtained with qir/tt = 0.93 and (in GeV): 
m = 0.007, A ir = 0.24, A uv = 0.905 [IJ. The Bethe- 
Salpeter amplitudes are canonically normalised; n-n is the in- 
pion condensate [5l| - [53| ; and f n , P are the mesons' leptonic 
decay constants. Empirical values are K n « (0.22 GeV) 3 and 
[H f n = 0.092 GeV, f p = 0.153 GeV. All dimensioned quan- 
tities are listed in GeV. 





F„ 


E P 


1/3 


m p fir 


fp 


3.639 


0.481 


1.531 


0.368 0.243 0.140 


0.929 0.101 


0.129 



forbidden by the interaction, then the rainbow-ladder 
vector-meson Bethe-Salpeter amplitude takes the form 



(17) 



where P^ 



0j 7 J + la = 7m- We assume isospin 



symmetry throughout and hence do not explicitly include 
the Pauli isospin matrices^ 

Values of some meson-related quantities, of relevance 
herein and computed using the contact-interaction, are 
reported in Table fl] We quote pion properties in order 
to provide a broader picture: the pion's Bethe-Salpeter 
amplitude is 



r*(P) = 76 



iK(P) + ^7 • PF^P) 



(18) 



III. p-MESON FORM FACTORS 

At this point we can proceed to computation of the 
form factors. In order to ensure a symmetry-preserving 
treatment, one must calculate the vertex in Eq. ([3]) at 
the same level of approximation as used for the dressed- 
quark propagator and meson Bcthe-Salpctcr amplitude; 
i.e., the generalised impulse approximation: 

T a ^(p,p') = Tl^{p,p>) + r d a ^(p,p') , (19) 
r^(p,p')=2 J 0^Tr CD l [ iT p /(k;-p')S(k ++ ) 
xzr/(fc_+, fc ++ )5(fc_+)zr« (k - 9 /2;p)S(fc__)}, (20) 



wherein the trace is over colour and spinor indices and 
k af3 = k + aq/2 + Pp/2. We illustrate Eq. @® in Fig.HJ 
In evaluating Eq. (fT9)) we write: 



S f = S + 5^fL Sf,f = u,d, 



(21) 



1 Note, too, that we use a Euclidean metric: {7^,7,/} = 28 ^ 
ll = 1 P ; 75 = 74717273, tr[757 fl 7, / 7 p 7 CT ] = -4e fl „ pCT ; cr M „ 
(«/2)[7m,7^]; a ■ b = / l^i ajbi; and P M timelike =>• P 2 < 0. 




k. 



FIG. 1. Impulse approximation to the p-7 vertex, Eq. (|20|l : 
solid lines - dressed-quark propagators; and shaded circles, 
clockwise from top - Bethe-Salpeter vertex for quark-photon 
coupling, and Bethe-Salpeter amplitudes for the p + -meson. 



where S is given in Eq. ((9j) , with the dressed-mass ob- 
tained from Eq. (fTTj|) . and the broken-CP corrections 
8gp Sf are detailed below; and the p-amplitude 



(22) 



with E p {P) explained in connection with Eq. (JTTJ) and 

p CP 

the broken-CP corrections Ta explained below. Our 
computed values for the dressed-quark mass, M, and E p 
are listed in Table Q] 

The remaining element in Eq. (|19p is the dressed- 
quark-photon vertex. We are only interested in the 
q 2 = values of the form factors and hence may use 

er M (pi,p 2 ) = eQ7 fl +i<D'y 5 (T )J , v (p2-pi) v (23) 
=: e diag[e u r^ ( Pl ,p 2 ), -e^ fa , p 2 )] , (24) 

where e is the positron charge, Q, = diag[e u = 2/3, — e<j = 
1/3] and T> = diag[d„, —dd], with df the EDM of a cur- 
rent quark with flavour /. N.B. The second term in 
Eq. (|23|) describes the explicit current-quark EDM inter- 
action in Eq. @. In Sec. lIVI we show that the other terms 
in Eq. ^) generate additional contributions that interfere 
with this explicit term. 

Note that both structures in the vertex, Eq. (|23|) , are 
in general multiplied by momentum-dependent scalar 
functions. Naturally, the vector Ward-Takahashi iden- 
tity ensures that the coefficient of the Qff^ term is "1" 
at q 2 = 0. In connection with the tensor term, one 
knows from Ref. [2(| that a tensor vertex is not dressed 
in the rainbow-ladder treatment of the contact inter- 
action. However, with a more sophisticated interac- 
tion, strong interaction dressing of the 750^ part of the 
quark-photon vertex might be significant, given that the 
dressed-quark-photon vertex certainly possesses a large 
dressed-quark anomalous magnetic moment term owing 
to dynamical chiral symmetry breaking j55j . At q 2 — 0, 
this could enhance the strength of the (D term by as much 
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as a factor of ten. If so, then sensitivity to current-quark 
EDMs is greatly magnified. It is worth bearing this in 
mind. 

Working with Eq. ([3]), it is sufficient herein to employ 
three projection operators: 



(25a) 



6p 2 



p3 



1 



with p' = p 



2iq 
q, P 



2 (p Ia>(p) £ c t >P>v.oqa ( Pj>l3(p') 

p' , for then 
1 



£(0) = lim 

q 2 ^o 12m 2 a 



1 

lim - P r 

9 2 ->o 4 

lim Pf 



fW(0) 
2>(0) 



(25b) 
(25c) 

(26a) 

(26b) 
(26c) 



and Hp = e/[2m p ], d p = 2)(0) e/[2m p ]. So long as a 

global-symmctry-prcscrving regularisation scheme is im- 
plemented, £(0) = 1; the value of !M(0) is then a pre- 
diction, which can both be compared with that produced 
by other authors and serve as a benchmark for our pre- 
diction of £>(0). 

At this point one has sufficient information to calcu- 
late the p-meson's magnetic moment. We simplify the 
denominator in Eq. (|19[) via a Feynman parametrisation: 

{k 2 ++ + M 2 ) ~* (k 2 _ + + M 2 ) 1 {k 2 __ + M 2 ) 1 



= 2 



dx dy 



o Jo 
1 



k 2 + M 



[p 2 -2(l-2x-2y)p-q + q 2 ] 



-(1 -2y)q-k + (l- 2x)p ■ k 



(27) 



This appears as part of an expression that is integrated 
over four-dimensional fc-space. The expression is simpli- 
fied by a shift in integration variables, which exposes a 
denominator of the form l/[fc 2 + M 2 ] 3 , with 



M 



M 2 +x(x - l)m 2 + y(l 



y)< 



(28) 



One thereby arrives at a compound expression that in- 
volves one-dimensional integrals of the form in Eq. (fTOjl . 
which we regularise via Eq. (JTTJ) and generalisations 
thereof; viz., 



ds 
ds 
ds 



S + UJ] 

s 



d_ 

do: 
\£_ 
2duj 2 



-i-.C™(u)=:C?(co), (29a) 
<?>)=: (29b) 
Ci"(w) -u7%{w), (29c) 
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FIG. 2. Evolution of p-meson magnetic moment with current- 
quark mass. 77i = 170 MeV corresponds to the mass of the 
s-quark in our treatment of the contact interaction so 
the difference between t M P {Q) and iW^(0) is just 1%. 



etc. Details for this component of our computation may 
be found in Ref. j20| and pursuing it to completion one 
obtains the magnetic moment listed in Table [TTT 

We depict the evolution of with current-quark 

mass in Fig.fSJ £4f(0) is almost independent of m. 
This outcome matches that obtained in Ref. (Hf using a 
renormalisation-group-improved one-gluon exchange ker- 
nel and hence a momentum-dependent dressed-quark 
mass- function of the type possessed by QCD (60l - [63j . The 
behaviour in Fig. [2] will serve to benchmark that of the 
p-meson's EDM. 



IV. p-MESON EDM: FORMULAE 

We now turn to computation of the effect of the inter- 
action terms in Eq. ^ on the p-meson. There are three 
types of contribution, which arise separately through 
modification of: (1) the quark-photon vertex, Eq. 



TABLE II. Magnetic moment of the p-meson calculated using 
our framework; and a comparison with other computations. 
Legend: RL RGI-improved, treatment of a renormalisation- 
group-improved one-gluon exchange kernel in rainbow-ladder 
truncation; EF parametrisation, entire function parametrisa- 
tion of solutions to the gap and Bethe-Salpeter equations; and 
LF CQM, light-front constituent-quark model. The results 
are listed in units of e/[2m p ]. 



This work and Ref. [20] 


2.11 


DSE: RL RGI-improved [34] 


2.01 


DSE: EF parametrisation [56] 


2.69 


LF CQM [57] 


2.14 


LF CQM [58] 


1.92 


Sum Rules [59] 


1.8 ±0.3 


point particle 


2 
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FIG. 3. Top - Correction to the quark-photon vertex gener- 
ated by the four-fermion operator in Eq. I|30p . The unmodified 
quark-photon vertex is the left dot, whereas the right dot lo- 
cates insertion of Lg. If the internal line represents a circulat- 
ing d-quark then, owing to the Is insertion, the external lines 
are it-quarks, and vice versa. Middle - Analogous correction 
to the p-meson Bethe-Salpeter amplitude. The unmodified 
amplitude is the left dot, whereas the right dot locates inser- 
tion of L%. The lower internal line is an incoming d-quark 
and the upper external line is an outgoing w-quark. Bottom - 
Xe-correction to the dressed-quark propagator, with the dot 
locating the operator insertion. If the outer line is a u-quark, 
then the internal line is a d-quark; and vice versa. 



(2) the p-meson Bethe-Salpeter amplitude, Eq. (|2"2")l ; and 

(3) the dressed-quark propagator, Eq. (|2"Tj) . 



Wick contractions produces the following result: 



'A 2 e„ 



(^[/, 12 + iV c /3], (31a) 



I 12 = -P R S(£ + q)^S(£)P R 

+P L S{£ + q)^S(£)P L , 
I* = P L tr{S(e + q) lfi S(£)P L } 



(31b) 



-P R ti-{S(£ + q)^S(£)P R }, (31c) 

where Pr,l = (1/2) (1 ±75). These right- and left-handed 
projection operators satisfy P R + Pl = Id- 
Further simplification of the integrand reveals 
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M 
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so that one may subsequently obtain 
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(32b) 
(32c) 



4 M 



(27T) 



(2tt) 4 " 
d A £ - 



iM 
16-k 2 



dx d (u q ) , 







iM 
iM 



"9m75^-2 / dxxd (U q ) 

^ Jo 



(33a) 
(33b) 



(2tt) 



4 M 



dx (1 - 2x) d K) , (33c) 



where 



Eq. (|31a[) becomes 



x(l — x)q 2 + M 2 . Combining the terms, 



A. Four-fermion interaction 

We begin with the dimension-six operator, which can 
be written explicitly as 



Le = i^ [u a d a S>^u h + u a lb d a d 



-d a d a u b lb u b - d a lb d a u b u b ] , (30) 



with summation over the repeated colour indices. This 
operator generates all three types of modification. 



%.e d M 



dx d (0Jq) 



A 2 e u 16tt 2 Jo 
x [(1 + 2JV C )(1 - 2x)q^ + Ur^qvYfr (34) 



g 2 =o X. e d M - iu 2 
= A^^16^ Cl(M) ^ 75 ' 



(35) 



In the other case, with a u-quark circulating in the 
loop, one obtains 



,7xff q 2 =o %. e u M -h 



Ta e * = " ^ — ^C^{M 2 )ia^q vl ^. (36) 



A 2 e d 16tt 2 

Plainly, the net correction to the quark-photon vertex 
can now be cast in the form of the second term in Eq. (|23|) 
and hence is readily expressed in 2)(0). 



1. Lq - quark-photon vertex 

This contribution is depicted in the top panel of Fig. [3] 
Consider first the case of rf-quarks circulating in the loop, 
then straightforward but careful analysis of the induced 



2. Le - Bethe-Salpeter amplitude 

This correction is depicted in the middle panel of Fig. [3] 
Each of the four terms in Eq. (|3U|) generates a distinct 
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contribution. That from the first and second are: 
T p a 4 = -i^N c E p P R 



xtr 



S(£)P R S(£ + P)~fZ, (37a) 



(27T)' 



■S(£ + P)^S(£)P R . (37b) 



The third and fourth terms are identical, up to sign- 
change and the replacement P R —> Pt,; and hence 

^ = i§ E P J [^ 2T + ^4 3T ] , (38) 

where the superscript "T" indicates that 7J is here used 
in the expressions for I 12 , I 3 . 

Now, using the formulae of Sec. lIV ATI one arrives at 



■ K ME p 
''A 2 16tt 2 



j5<r a vPv / dxCT(u P ), (39) 



where wp — x(l — x)P 2 + M 2 , P 2 = —m 2 . This is one of 
the additive corrections to the Bethe-Salpeter amplitude 
anticipated in Eq. 




■9/2 



P + 9/2 

FIG. 4. Correction to the quark-photon vertex generated by 
the quark chromo-EDM operator in Eq. © : the incoming and 
outgoing quark lines have the same flavour, /. The dot in the 
left two diagrams locates insertion of Lcedm, whilst that in 
the rightmost diagram indicates the second term in Eq. (|23l) ; 
i.e., the explicit quark EDM. 



B. Quark chromo-EDM 

The term in the middle line of Eq. ([2j also generates 
all three types of modification described in the opening 
lines of this Section. Notably, owing to dynamical chi- 
ral symmetry breaking, the dressed-quark-gluon coupling 
possesses a chromomagnctic moment term that, at in- 
frared momenta, is two orders-of-magnitude larger than 
the perturbative estimate [55|. One may reasonably ex- 
pect similar strong-interaction dressing of a light-quark's 
chromo-EDM interaction with a gluon, in which case sen- 
sitivity to the current-quark's chromo-EDM is very much 
enhanced. 



3. Lq - quark propagator 



The final modification arising from the dimension-six 
operator is that depicted in the bottom panel of Fig. [3] 
So long as the correction is small, it modifies the dressed- 
quark propagator as follows: 

S{k) -> S(k) + Sj^Sik) = S(k) + S(k)iT SLe S(k) , (40) 



1. La 



quark-photon vertex 



This contribution is depicted in Fig.|3] After a lengthy 
analysis, in which we represent the exchanged gluon via 
Eq. © , the sum of the two leftmost diagrams produces 



where, once again, each of the four terms in Eq. (|30j) con- 
tributes. Their sum is 

r S£e = 5/t04 [PrS(£)Pr PlS(£)P l 



(27T) 

-N c P R tr{S(£)P R } - N C P L tr{S{e)P L }] . (41) 



Now 



PrS(£)P r - P L S(l)P L 



AI 



P + AP 



75 



; [P R tv{S(£)P R } - P L tr{S(£)P L }] , (42) 



p7(ff) 



p bin m G J 

x ^2q a 07^75 - 6i[3(x - 1/2)^ - p u ] 75 

m G Jo I 
-6[(x - l/2)ojg + 2x(l - x)q-p] g M 7 5 

+M[({x - l/2)q+p) ■ 7\q a a afl j 5 
+Mq a o- a ^ 5 [((x - 1/2)9 + P) ■ 7] 1 , (44) 



so that with little additional algebra one arrives at 



5 u S(k) 



k 2 + AP 



(i+2Ay 5i^ iu(M2)75 - (43) 



where, again, d/ is the chromo-EDM of a quark with 
flavour /. 

As we are interested solely in the EDM, we may con- 



p-P/2 _ p - P/2 




e + p/2 



p-i 



V + P/2 




e + p/2 



p + P/2 



FIG. 5. Correction to the p-meson Bethe-Salpeter amplitude 
generated by the quark chromo-EDM operator in Eq. ([2]): the 
incoming line is a d-quark and the outgoing line is a tt-quark. 
In each case the dot locates insertion of Lcedm- 



+ 



FIG. 6. Correction to the dressed-quark propagator generated 
by the quark chromo-EDM operator in Eq. @. In each image 
the dot locates insertion of Lcedm- 



obtained from Eq. (|46[) via the interchange d u f-> dd, and 

p — > —p, P -> —P. 



3. Lce dm - quark propagator 



sider q 2 = 0, at which value the result simplifies greatly: 

p 2 ( 9) = l!W \c™(M 2 ) - C?(M 2 )] 
p 3i7r m G 



X [75 0> + 3ip M 75 
1 dfa m — i 



37T TOg 



d (M ) 



M{7-p,7 5 cr /ja }g Q 



+2p • gg M 7 5 + 6M 2 p^ 5 



(45) 



Plainly, the net correction to the quark-photon vertex 
from these two diagrams can now be cast in the form of 
the second term in Eq. (|23[) . which, in fact, is precisely 
the rightmost diagram in Fig. 2] because q =pi — p\. 



The last modification generated by the chromo-EDM 
term in Eq. is that to the quark propagator, Fig. [5] 
The self-energy insertion is readily evaluated: 



r*oo = d f - ^ © iu (M 2 ) 7s 



where 



(47) 



s 2 r T " 2 

£> m (w) = / ds j~^ I *r — exp(-rw), (48) 



so that, with / = u, d, 

W) = ^^/^^ 2 h5. (49) 



2. Lcedm - Bethe-Salpeter amplitude 

This correction is expressed in Fig.[SJ Owing to sim- 
ilarity between the -L e ff-uncorrected p-meson amplitude 
and quark-photon vertex, the results can be read from 
those in Sec. lIVBT] viz., with d± = d u ± dd, 



pp(s) 



1 «IE 



E p [ dx[C u {up)-C?{u P j\ 
Jo 

I [(d+ - 3{x - l/2)d-)P p - d-pp] *^ 5 <P f 



xhd+[u P -2M 2 ] Ppl ^ 



T 



-J_([cj P - 2M 2 ][{x - l/2)Pp+pp\) ilso- p ^ 
+Md d [((x - l/2)P + p)-'y]Ppo-p^ 5 2' f 



T 
f.ia 



+Md u Ppap li75 ^ a [((x - l/2)P + p) • 7] I • (46) 

In computing the vertex in Eq. (|19p one must employ 
Fig. [5] and also its charge conjugate, the form of which is 



C. 6>-term 

Owing to a connection between the Higgs mechanism 
for generating current-quark masses in the SM and CP 
violation in the weak interaction, the effect of the #-term 
can completely be expressed through a Ua(1) rotation 
of the current-quark mass-matrix. We consider the s- 
quark to be massive and m u — m,d, in which case the 
effect of the first term in Eq. is expressed simply in a 
modification of the dressed-quark propagator: 



S{k) 



17 • k + M + |m07 5 



? small 



S(k) 



1 



k 2 + M 2 2 



rati 75 



(50) 



(51) 



1. Dressed-quark anomalous chromomagnetic moment 

In our global-symmctry-prcscrving rainbow-ladder 
treatment of the contact interaction, the general form 
of the p-meson's Bethe-Salpeter amplitude is given in 
Eq. (fTT)) . The absence of a term a flv P u F p (P) is an arte- 
fact of the rainbow- ladder truncation: even using Eq. (JS]) , 
a Bethe-Salpeter amplitude with F p (P) ^ is obtained 
in any symmetry-preserving truncation that goes beyond 
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this leading order [3C§ . One material consequence of this 
omission is complete cancellation of all terms at leading- 
order in 9, so that the ^-term's contribution to the p- 
meson's EDM is anomalously suppressed in rainbow- 
ladder truncation. This defect may be ameliorated by 
acknowledging that the dressed-quark-gluon vertex pos- 
sesses an anomalous chromomagnetic moment coupling 
which is enhanced by dynamical chiral symmetry break- 
ing [55[. We therefore include an effect generated by 



-pacm 



(PuPf) 



2AI 



(52) 



where [6J ^ acm - (-1/4). 

In order to explicate the effect we find it convenient to 
first express collectively the corrections to the dressed- 
quark propagator computed above; viz., from Eqs. (|43l) . 

m, dsn), 



S(k) 

Xj^ 



Sik) - 175 



A 



-(1 + 2JV C ) 



k 2 + AP 
K M 



A 2 16tt 2 



A 



(9) 



Xs 



~ 8 OIR 

7r 



-df - ^ © iu (M 2 ) , 



-m6 . 



(53) 

C iu {M 2 ), (54) 
(55) 

(56) 



Our corrections are now obtained via the diagrams in 
Fig-Ill except that here the dots represent Eq. (|52|). and 
one simultaneously adds the correction to one and then 
the other propagator. In this way, careful but straight- 
forward computation yields 



r 



A.acm 



~, \ 1 ,,acm \2 , ,ac 



2iitm 2 ^ 



2M 



x / dx [C iu (up) - C?{uj p )} 7m75 
Jo 

x|3 M i cm 7-(P+(^-l/2)P) 

x[(l -x)A 1 -xX 2 ]P,, 
+i [(1 - x)X 1 + xX 2 )} [rt cm l»P a a ali 

-X-M[^ m (p + (x-l/2)P)fi^p 

+3i^(p + (x - 1/2)P)„] | 75 , (57) 

where pJ± m = /i acm ± /i| cm , and {X i , i = 1, 2} represents 
the quark propagator correction on each leg with A ± = 
A X ±A 2 . 

One can now adapt the general expression in Eq. (|57[) 
to the particular cases of relevance herein. The first is 
the p-mcson Bcthc-Salpcter amplitude. Capitalising on 



isospin symmetry, which entails /i acm = H d cm =: /i acm , 
one finds 

acm\— pi 



X7M^a75 

am M° 
3iTrm 2 G 2M 



E p / dxC[ u {ujp)\i[{l~x)X 1 



+xX 2 ] ( lfi P v o va - P v o vc flfi) (p+(x- l/2)P){3 
-2X-M^ a (p+(x - l/2)P) tf <7^ I76 , (58) 



where "A" is constructed from the correction specified in 
one of Eqs. ([M]) - 1(550 ■ 

The other case is the quark-photon vertex, for which 
the correction is found with A 1 = A 2 = A, since the quark 
flavours are identical, and we need only consider q 2 = 0: 

n ..acm\ _ 

r ^=^^ r ^(M 2 M,. P ,a, aqa ]. (59) 



V. p-MESON EDM: RESULTS 

A. Analysis without Peccei-Quinn symmetry 

In order to obtain a result for the p-meson's EDM, d p , 
it remains only to sum the various contributions derived 
in Sec. lIVI as they contribute to Eq. (fT9|) . evaluated with 
the parameter values in Table |TJ 

d p = -2.88 x 10" 3 n &CTa eQ/s 
+0.785 {d u - d d ) 
+ (1.352 + 0.775 n acm )e{d u - d d ) 
-(0.091 - 2.396 n acm )e{d u + d d ) 

-e^y (2.696 - 6.798 /i acm ) x 10~ 3 . (60) 

In this formula, df,df carry a dimension of inverse- mass 
and s = 1 GeV. 

A nugatory transformation allows one to rewrite 
Eq. (|60[) in terms of dimensionless electric and chromo- 
electric quark dipole moments; viz., 



d n = 



-2.88 x 10" 3 pT m ee/s 



A 2 



0.785 (D u - D d ) 



-(1.352 + 0.775 n acm )e(D u - D d ) 
-(0.091 - 2.396 /Lt acm )e(D„ + D d ) 

-(1.096 - 2.763 ^ acm ) x 10~ 5 eXi 



, (61) 



where Vh = 246 GeV is the cube-root of the phcnomeno- 
logical Higgs vacuum expectation value. In a class of 
models that includes, e.g., Ref. [65[, one finds 



D I 



m f 
vh 



2 x 10" 



(62) 



10 



o.90r 




0.75L , , , J 

50 100 150 200 

m (MeV) 

FIG. 7. Evolution of the quark-EDM component of the p- 
meson's EDM with current-quark mass, assuming d_ is inde- 
pendent of m. m — 170 MeV corresponds to the mass of the 
s-quark in our treatment of the contact interaction [z|, so 
the difference between (P p and is 10%. 



TABLE III. Contributions to the p-meson EDM associ- 
ated with a quark chromoelectric dipole moment, with 
d\ — e(d u =p dd). Row 1: quark- photon vertex correction, 
Sec. lIV B~T1 Row 2: p-meson Bethe-Salpeter amplitude correc- 
tion, Sec. lIV B~2l Row 3: dressed-quark propagator correction, 
Sec. lIV B~3l Row 4: anomalous chromomagnetic moment con- 
tributions, Sec. lIV C~I1 Row 5: sum of preceding four rows; 
Row 6: Row 5 evaluated with p acm = -1/4; and Row 7: 
sum rules result from Ref. [l5j |. evaluated here with a heavy 
s-quark. 



qjq 




-0.066 dt 


- 0.199 d+ 


BSA 




-0.120 dt 


+ 0.108 


S(k) 




1.538 dt. 




acm (xp acm ) 




0.775 dt 


+ 2.396 d+ 


our CEDM 


(1.35- 


-0.78p acm )dl - 


- (0.09 - 2.40p acm )d+ 


total 




1.16 dt 


- 0.69 d+ 


sum rules [151 




-0.13 dt 





a result which may be used to inform expectations about 
the "natural" magnitude of the terms in Eqs. ([60]) . (|61[) . 

There are four distinct types of contribution to d p in 
Eq. (|60[) . The first is associated with the #-term; and it is 
notable that this contribution vanishes in the absence of 
a drcssed-quark anomalous magnetic moment, a feature 
which emphasises the connection between topology and 
dynamical chiral symmetry breaking that is highlighted, 
e.g., in Eq. (21) of Ref. j4Q|. Our result may directly be 
compared with that obtained in a sum rules analysis; viz., 



herein 
Ref. [15] 



-2.9 x 10~ 3 A* acm e6» 
4.4 x 10~ 3 e<9. 



0.7 x 10- 3 e6> 



(63) 



The second contribution owes to an explicit dressed- 
quark EDM. It has been computed via a number of meth- 
ods, so that a comparison with our results is readily com- 
piled: 

herein DSE [13] BM [13] nrQM [13] sum rules [15] 
0.79 0.72 0.83 1.00 0.51, 

(64) 

where each entry is multiplied by d- = (d u — dd)] and 
DSE [l3j summarises results obtained from momentum- 
dependent DSE input, BM [l]| reports a bag-model re- 
sult, and nrQM [131 ] is the non-relativistic constituent- 
quark value. Wc depict the current-quark mass depen- 
dence of this contribution in Fig. [7] It is notable that the 
magnitude of these results matches an existing DSE esti- 
mate of the analogous contribution to the neutron's EDM 
|66j ] . Moreover, based on Ref. [g7j , a perturbative analysis 
would yield 2m p d p ' crt = 2m d„ , where m is the current- 
quark mass. With the parameter values employed herein, 
this is dP p CTt = 0.014 cL, which is just ~ 2% of the order- 
of-magnitude specified by the values in Eq. (|64p . 

The third contribution to d p is generated by the 
quark's chromoelectric dipole moment. Its subcompo- 



nents are detailed in Table IIIII The net result is com- 
parable in magnitude and sign with that produced by 
the quark EDM, Eq. (|64p . In comparison with a sum 
rules computation [151 ] . however, our result is an order of 
magnitude larger, has the opposite sign and contains a 
sizeable d+-term. At least the first two of these marked 
discrepancies are insensitive to reasonable variations in 
/i acm . It is worth emphasising here that our calculation 
has no other variable parameters: the two specifying our 
model, listed in Table U were fixed in prior studies of 
an array of meson and baryon observables [l7l - l22T ] . This 
mismatch will receive further attention in future work. 

The four-fcrmion interaction is responsible for the final 
contribution in Eq. (|60|) . Its subcomponents are detailed 
in Table IIVI As ours is the first estimate of the contri- 
bution from a dimension-six operator to the p-meson's 
EDM, there is no ready substantial comparison. On the 
other hand, the result in Table |IV] is quickly seen to be 
"natural" in size. The dimension-six operator is associ- 
ated with a coupling ^C/A 2 , which has mass-dimension 
"—2" . In order to obtain a quantity with mass-dimension 
"—1", this coupling must be multiplied by another en- 
ergy scale. We are interested in an hadronic EDM, so 
that scale should be typical of hadron physics; e.g., the 
dressed-quark mass "M". Finally, a loop correction is 
required for the generation of an EDM, and loops are 
characterised by a factor l/(167r 2 ). Putting these quan- 
tities together yields an expectation based on naive di- 
mensional analysis; viz., 



jD=6 



1 M Vh-K 



I671- 2 vh A 2 



1 x 10 



- 5 e ggjC 
A 2 



(65) 



in agreement with the magnitude of the final row in Ta- 
ble |lVl Comparison with Eq. ([62]). furthermore, indicates 
that in our computation the quark-EDM and dimension- 
six contributions are naturally related via 



(66) 
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TABLE IV. Contributions to the p-meson EDM associated 
with the dimension-six operator in Eq. Each row should 
be multiplied by evn'K/ A 2 . Row 1: quark-photon vertex cor- 
rection, Sec. lIV A~D Row 2: p-meson Bethe-Salpeter ampli- 
tude correction, Sec. lIV A~2l Row 3: dressed-quark propaga- 
tor correction, Sec. lIV A~3l Row 4: anomalous chromomag- 
netic moment contributions, Sec. lIV C~TI and Row 5: sum of 
preceding four rows. 



qjq 


-1.005 x 10" 5 


BSA 


-9.114 x 10~ 7 


S(k) 





acm (x^ acm ) 


2.763 x 10 _5 p acm 


our D = 6 total 


-(1.096 - 2.763 p acm ) x 10 -5 


B. 


Peccei Quinn Symmetry 



The leading term in Eq. (pTTj) is that associated with 9. 
Arising from a dimension-four operator, this contribu- 
tion is not suppressed by a large beyond-SM mass-scale. 
One may furthermore expect that, absent any symme- 
try to prevent it, a typical non-SM for CP- violation will 
produce large corrections to 9. In order to reconcile this 
with the remarkably small upper-bound on 9 placed by 
the neutron's EDM, one must accept that the initial value 
of 9 is very finely tuned. There is nothing to prevent this 
from being simply an accident of Nature. However, some 
view that possibility as aesthetically displeasing and pre- 
fer to introduce a new dynamical degree of freedom, the 
axion, a pseudo-Goldstone boson, whose role is to cancel 
the effect of 9 [(38| . It is notable that there is currently no 
empirical evidence in favour of the axion's existence and 
the remaining domain of parameter space is small j69| . 

Notwithstanding this, in the context of EDM estimates 
it is customary to expose the possible effect of axion 
physics on the results in Eq. or (|6 1 [) . Here there 
is a complication. If one considers an extension of the 
SM with a collection of CP-odd operators that may mix 
with the 9-teim, then the effective potential describing 
axion physics at the hadronic scale can plausibly acquire 
terms that shift its minimum to a nonzero value of the 
effective ^-parameter, induced 0- The quark chromo- 
electric dipole moment interaction is one such operator. 
In its case, within a sum rules calculation fTo| . the net 
effect of this mixing is elimination of 9 in favour of a 
modest enhancement in magnitude of the coefficients of 
d± in Eq. (|60[) . with no change in sign. 

The implications for our study are plain. Allowing 
an axion-like mechanism to play a role, then 9 disap- 
pears from Eqs. (|60[) and (|61|) . and any measurement of 
an hadron EDM, here that of the p-mcson, places a little 
more stringent constraint on d± in particular but also on 
d± and 

This is, perhaps, particularly relevant to since the 
high-scale physics that generates this operator will typi- 
cally also produce a complex phase for the quark masses. 



Within the low-energy effective theory of Eq. ([2]) , this 
phase will arise from one-loop contributions to the quark 
propagator containing one insertion of the CP-violating 
four-quark operator and the quark Yukawa interaction. 
Consequently, constraints on 9 imply a bound on On 
the other hand, with the elimination of 9 via an axion 
effective potential, the term modulated by %^ is exposed 
to independent constraint Q. Computing the contribu- 
tion of the four-quark CP-violating operator to the ax- 
ion potential, determining the resulting dependence of 
^induced on and deriving the expression corresponding 
to Eq. (|5T|) will be the subject of future work. 



VI. EPILOGUE 

Using the leading-order in a global-symmetry- 
preserving truncation of QCD's Dyson-Schwinger equa- 
tions, we computed the electric dipole moment of the 
p-meson, d p , that is generated by the leading dimension- 
four and -five CP-violating operators and an example of 
a dimension-six operator. We employed a momentum- 
independent form for the leading-order kernel in the 
gap- and Bethe-Salpeter equations. This is known to 
produce results for low-energy pseudoscalar- and vector- 
meson obscrvables that are indistinguishable from those 
obtained with the most sophisticated interactions avail- 
able when they are analysed using the same truncation. 
Since the dipole moment is a low-energy observable, our 
predictions should be similarly reliable, in which case 
the framework we employ and elucidate can usefully be 
adapted to the more challenging task of computing the 
neutron's EDM, d n . 

We find that the two dimension-five operators; namely, 
quark-EDM and -chromo-EDM, characterised by d q and 
d q , respectively, produce contributions to d p whose coef- 
ficients are of the same sign and within a factor of two 
in magnitude. This contrasts with an extant sum rules 
evaluation, in which the coefficients of the contributions 
have the opposite sign and differ by a factor of four in 
magnitude. Since all studies agree within a factor of two 
on the quark-EDM coefficient, the discrepancy resides 
with the chromo-EDM contribution. These differences 
invite further analysis and guarantee relevance to a DSE 
evaluation of the impact of d q on the neutron's EDM. 

Absent a mechanism that suppresses a 0-term in 
any beyond-Standard-Model action, the tight constraints 
on the magnitude of a contribution from this term to 
the neutron's EDM also apply to contributions from a 
dimension-six four-fcrmion operator to this or another 
hadron's EDM. Should such a mechanism exist, how- 
ever, we find that a dimension-six operator can match 
the quark-EDM and chromo-EDM in importance. 

Using the techniques described herein, calculation of 
the neutron's EDM is underway. 
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